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The generalized Swiss - cheese model, consisting of a Lemaitre - Tolman (inhomogeneous dust) 
region matched, by way of a comoving boundary surface, onto a Robertson- Walker background of 
homogeneous dust, has become a standard construction in modern cosmology. Here we ask if this 
construction can be made more realistic by introducing some evolution of the boundary surface. The 
answer we find is no. To maintain a boundary surface using the Darmois - Israel junction conditions, 
as opposed to the introduction of a surface layer, the boundary must remain exactly comoving. The 
options are to drop the assumption of dust or allow the development of surface layers. Either option 
fundamentally changes the original construction. 
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I. INTRODUCTION 

In [l| we generalized the Swiss-cheese cosmologies so 
as to include nonzero linear momenta of the associated 
boundary surfaces. We found that the final effective grav- 
itational mass and size of the evolving inhomogeneities 
depends on their linear momenta but these properties are 
essentially unaffected by the details of the background 
model. In this paper we explicitly specify the inhomo- 
geneity and take it to be Lemaitre - Tolman (L - T) 
inhomogeneous dust. With a comoving boundary sur- 
face, this construction is a standard feature in modern 
cosmology Since an exactly comoving boundary 
surface is clearly a mathematical idealization, the pur- 
pose of this paper is to examine possible evolution of the 
boundary surface. We find that if we rigorously maintain 
the Darmois - Israel junction conditions, no evolution of 
the boundary surface is possible. Hints that this stan- 
dard construction is a mathematical artifact are reviewed 
in Evolution away from dust, or the development of 
surface layers (which we briefly discuss) fundamentally 
changes the original construction. 



II. LEMAITRE-TOLMAN SOLUTIONS IN 
LAGRANGIAN FORM 

We find it convenient to use the somewhat less common 
Lagrangian coordinates to describe the L - T solution and 
we begin by summarizing some fundamental properties 
of the solution in these coordinates. 



A. Basic calculations 

Starting with Einstein's equations @ 

G a p + ^g a p = S>TrT af 3, (1) 
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where A is the cosmological constant, we seek solutions 
for dust, 



T a /3 = P UaU/3 



(2) 



where the u a are tangent to the generators of the geodesic 
flow and we consider only positive definite energy densi- 
ties p > 0. We use comoving synchronous coordinates so 
that 

ds 2 = e a(m ^dm 2 + R 2 (m, t)dtt 2 - dt 2 (3) 

where dVl 2 is the metric of a unit two-sphere, which we 
write in the usual form d9 2 + sin 2 (9)d(f> 2 , and we assume 
the existence of an origin defined by R(0, t) = 0. The 
generators of the flow are then u a 
normals are n a = ±e~ a / 2 5^ so that 
and u a n a = 0. From 



<5" and the radial 

\ a u n = n a n n = 1 







we find 



(R) 



'\2 



1 + 2E(m) 

where E is an arbitrary function (> —1/2). Define 



(4) 



(5) 



(6) 



where 1Z is the Riemann tensor and so M. is the effective 
gravitational mass [l| . We obtain 



R 2 



2E 



2M 



With the gauge condition 



R d 

M(m,t) =m + A — 
6 



(7) 



(8) 



the information required to solve Einstein's equations is 
simply 0, the m derivative of © and equation ([8}. We 
arrive at 
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47rp(m, t) 



R 2 R' 



(9) 
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The implicit solution to §7$ is 

'* dx 



E. Conformal Flatness 



t-t B {m) = ± f 
Jo 



(10) 



J2E+±p + m 



whcrc Ib is a second arbitrary function. 



B. Invariants 

The L - T solutions have two independent invariants 
derivable from the Ricmann tensor without differentia- 
tion. These can be taken to be 



and 



K = 2(471-/9 + 2A) 



16 / 3m 



(12) 



where 1Z is the Ricci scalar and w is the first Weyl in- 
variant (C a p-ysC a ^ lS where C a p^s is the Weyl tensor). 



C. The Origin 



From ([9]) we have 



R ' 3 

lim — = lim 3R 2 R = - j- — r. 

m->0 m rn->0 47rp(0, t) 



and so from (TT21 

w(0,t) = 0. 

Further, it follows immediately from ([7|) that 



lim E = 0. 

m— >0 



(13) 



(14) 



(15) 



D. Singularities 

It is clear from © and (jlll) that scalar polynomial 
singularities occur for 



i^i? =0 



(16) 



where p, and therefore 1Z and u>, diverge. A "bang" (or 
"crunch") occurs for R = 0. Shell crossing singularities 
occur for R = 0. The conditions for their avoidance 
are well known for A = [|| and these conditions have 
recently been extended to the case 



Conformal flatness of the L - T solutions can be recog- 
nized by the following necessary and sufficient condition 
for conformal flatness [ijj: A = where 



A = m ttR 3 o. 

3 H 



From (HTJ) and © then 



R 



R 
3m 



so that 



R = a(t)m 



1/3 



( n ) Using (0 in Q with ||SJ) we find 



2E 

m 2/3 



a 2 A 



(17) 



(18) 



(19) 



(20) 



and so both sides of this relation have to be constants 
(the left being independent of t and the right being inde- 
pendent of m). The constancy of the term on the right 
is, of course, the Friedmann equation for dust. Finally, if 
a(t = c) = it follows from (fTU)) that ts(m) = c where c 
is a constant. In summary, the condition A = reduces 
© to the Robertson - Walker (R - W) metric. 



III. JUNCTION CONDITIONS 

The theory of hypersurfaces (S, here considered to be 
timelike) in spacetime is well established (see, for exam- 
ple, [ll|) but we think it important to cover a few details 
here. We follow the notation of [l2j. The Darmois [HI - 
Israel [l4| junction conditions for boundary surfaces re- 
quire the contin uity of the first and second fundamental 
forms accross E [l5| ■ These conditions are in general not 
equivalent to the conditions of Lichnerowicz [l6j which 
require the continuity of the metric and all first order par- 
tial derivatives of the metric across E in coordinates that 
traverse E (such coordinates being referred to as "admis- 
sible" ) . The only case when the two sets of conditions are 
equivalent is when Gaussian - normal coordinates (i~7| are 
used. For example, taking (m, 9, cj>, t) continuous through 
E, and using the Lichnerowicz conditions, we arrive at 
the constraint [p] = even for dm/dr = along E. This 
constraint is erroneous as we discuss below. Whereas one 
might argue that the coordinates (m, 0, <j>, t) are not "ad- 
missible" , here, to avoid any source of confusion, we rely 
only on the Darmois - Israel conditions. 

The basic model consists of randomly distributed non- 
intersecting timelike spherical boundary surfaces E in a 
R - W background V + . The local inhomogeneities V~ 
arc L - T solutions. Both V have metrics of the general 
form with V + being conformally flat. On E we write 



gijdx l dx 3 



(21) 
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where gij is the metric intrinsic to E (first fundamental 
form) and x 1 are the coordinates intrinsic to S. With 
spherical symmetry, without loss in generality we take 
the coordinates 9 and <fi (but not t) continuous through 
E on which we write the intrinsic metric as 



dsl 



R 2 (r)dn 2 - dr 2 



(22) 



where r is the proper time on E. Note that the proper 
time derivative d/dr is defined only along E. 

Now our use of the Lagrangian variable m requires 
special attention. A straight - forward calculation gives 



K 2 




(23) 



where Kij is the curvature extrinsic to E (second fun- 
damental form). For boundary surfaces, if follows from 
((2"3")) that [Ai] = 0. Throughout we consider A a uni- 
versal constant and set [A] =0. It now follows from © 
that 



0. 



(24) 



We write in"!^ as the unit radial 4 - orthogonal to E 
tangent to E. We simply drop 



and as the unit 4 

1 2 and write 



and 



-a/2* o e a/2*^ 

dr dr 



dm dt 
dr ' ' ' dr 



(25) 



(26) 



The vectors h a and u a are unit vectors in consequence 
of the relation 



dm 
~d~r~ 



-1 



(27) 



which follows from the fact that r is the proper time on 
E. We choose the global time-orientation dt/dr > 0. The 
use of the Lagrangian variable m makes the rigging of n a 
straightforward as m increases monotonically away from 
the origin. When E is comoving dm/dr = and we get 
fi a = n a and u a = u a , the normals and generators of the 
fluid flow. 



A. Boundary surfaces 

Let us assume that E is a boundary surface. Then, by 
assumption, the following well - known necessary (but 
not sufficient) conditions hold: 

[G a0 n a fi ] = [G a/3 u a n p ] = 0. (28) 

From Einstein's equations then we have 

[A] = [B] = (29) 



where 



and 



B = T a au a n . 



From (P), © and © we find 



if) 



and 



B = e a l 2 p— — . 

dr dr 



(30) 
(31) 

(32) 
(33) 



The standard assumption is to set E comoving so that 
dm/dr = and as a result the conditions (|29|) arc empty 
in the sense that A = B = 0. However, in this case, 
in addition to (|24[) . one very frequently encounters the 
additional constraint [p] = (which, as explained above, 
follows from the Lichnerowicz conditions) so that from 
((17)1 we have [A] = and so A = along E. To see 
that this additional constraint is spurious, one need only 
reflect upon classics: the Einstein - Straus model or Op- 
penhcimer - Snyder collapse! 



B. Non-comoving boundary surfaces: dm/dT\s 7^ 



From (El) and (US) we have 



" B 2 ' 






B 2 - A 2 _ 







From (jMJ) and [u a u a ] = then 

dm 
~dr~ 



so that from J32 



[p] = 0. 



(34) 



(35) 



(36) 



Condition ()36t is, for dm/dr^ 7^ 0, not an additional 
constraint, but rather a requirement of the Darmois - 
Israel junction conditions. In summary, with (|36p . for a 
boundary surface with dm/dr^ 7^ 0, and (| 1 7[) we have 
shown that 



[A] = 0. 



(37) 



This tells us that if E is a boundary surface and 
dm/dr |s ^ along E, then V~ is conformally flat at 
E. Now since dm/dr |s 7^ we can consider a region 
J in V~ defined by the non-vanishing ranges 5m and St 
bounded by E, to the future of E for dm/dr\^ > and to 
the past of E for dm/dr\s < 0. Since the unique evolu- 
tion of an R - W flow is the same R - W flow, the R - W 
solution in V~ is the same as that in V + . There cannot 
then exist non-comoving boundary surfaces in spherical 
dust that signify local inhomogencitics. 
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C. Surface layers 

The Israel formulation of thin shells (which we discuss 
only briefly here for completeness) involves the construc- 
tion of surface properties of S out of discontinuities in 
the extrinsic curvature, by way of the Lanczos equations 
[TT| fl2| [l4j . In particular, the surface energy density of 
£ is given by 



ina = —7 — "fijU l u 3 , 



(38) 



where 7^ = [Kij], 7 = 7^ and u % (= dx l /dr) is the three 
- tangent to £. The surface pressure (p) is defined by 



16np = 7 — jijjj^vP 



(39) 



These phenomcnological properties of £ are related to 
the evolution of £ via 

[n a u p V f3 u a ] = -7yuV = 8tt( P + a/2), (40) 

and to the enveloping 4 - geometries via 

u l V t a + (a + p)\7 iU l = [B] (41) 

and 



(a + p)n a v,P\7 pu a + pK\ = [A] 



(42) 



For the cases under consideration here we continue to use 
the Lagrangian variable m and continue to classify the 
motion via dm/dr^. 

For shells when we write down the standard equation 
for £, dm/dr^ = 0, we must keep in mind that [m] is 
not necessarily 0. By writing down dm/dT\s = we 
mean dm/dr]^. — dm/drlj = (and so [m] is constant 
along £). It then follows immediately form §3§ that £ is 
geodesic in V*. As a result, it follows from (|40|) that 



P- 



? = 



(43) 



At first glance, the possible cases are: p = a = 0, the 
boundary surface discussed above, and the thin shell with 
surface equation of state p — —a/2. However, from ([2]) 
and f4"2"j) we have 



pK$ = 



(44) 



and so either p = (again the boundary surface) or Kg 
changes sign across £. Since 



R 3 R 

?a/2 



(45) 



it follows that R , and from ([9]) also p, must change sign 
across £. We reject this possibility and conclude that £ 
cannot represent a comoving thin shell. 



Now for dm/dr|s ^ we have two possibilities: (i) 
dm/dr|^ = dm/dr^ ^ (and so [m] remains constant 
along £) and (ii) dm/dr]^ 7^ dm/dr]^ (with at least 
one 0). Clearly the notation (dm / dr\s) is inadequate 
for the second case. There is now a plethora of possi- 
bilities to consider for both cases, some of which have 
been examined in detail (see, for example, Q and refer- 
ences therein). Basically, we can choose some evolution 
for £ and see what surface properties (a and p) result, 
or choose these surface properties and see what evolu- 
tion of £ results. In both cases, the assumption that 
we have enveloping dust enters by way of the constraint 
equations (|4"Tj) and (|42p . To represent an inhomogeneous 
cosmological model, as opposed to a purely mathematical 
construction, physical input for a and p, or the evolution 
of £, is needed. In any event, we are now quite some 
ways away from the generalized Swiss - cheese model. 



IV. CONCLUSION 

We have examined the generalized Swiss - cheese 
model, consisting of a Lemaitre - Tolman (inhomoge- 
neous dust) region matched, by way of a boundary sur- 
face, onto a Robertson- Walker background of homoge- 
neous dust. A standard construction in modern cosmol- 
ogy is to assume that the boundary surface is exactly 
comoving. Here we asked if this construction can be 
made more realistic by introducing some evolution of the 
boundary surface. The answer we have found is no. We 
have found that to maintain a boundary surface using 
the Darmois - Israel junction conditions, as opposed to 
the introduction of a surface layer, the boundary must 
remain exactly comoving. We conclude that this stan- 
dard construction is a mathematical artifact, not an ac- 
ceptable physical model of a cosmological inhomogencity. 
The options are to drop the assumption of dust or allow 
the development of surface layers. Either option funda- 
mentally changes the original construction. 

In closing we should note that the situation we have 
considered, though common, is very idealized in the sense 
that the background has been taken to be exactly Robert- 
son - Walker. If this is not the case, £ can indeed evolve. 
Further discussion of various scenarios which have been 
studied can be found in H and [a. 
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